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Abstract 

We study a simple model of periodic contraction and extension of large intruders in a granular 
bed to understand the mechanism for swimming in an otherwise solid media. Using an event-driven 
simulation, we find optimal conditions that idealized swimmers must use to critically fluidize a sand 
bed so that it is rigid enough to support a load when needed, but fluid enough to permit motion 
with minimal resistance. Swimmers - or other intruders - that agitate the bed too rapidly produce 
large voids that prevent traction from being achieved, while swimmers that move too slowly cannot 
travel before the bed re-solidifies around them i.e., the swimmers locally probe the fundamental 
time-scale in a granular packing. 
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Granular materials are able to act either as a solid, supporting a load like sand on a beach, 
or as a fluid, e.g. flowing in an hourglass. Some lizards are known to make use of the balance 
between solid-like and fluid-like properties to swim through sand beds (Fig. 1). Since they 
cannot survive on the surface of their desert habitats during the day, they have developed 
the ability to submerge up to 10 cm beneath the sand to lower temperature regions, and 
even to travel within the sand bed. These reptiles exhibit a unique ability to manipulate the 
properties of granular materials, and their study stands to provide insight into the behavior 
of intruders in agitated granular beds more broadly, a topic of ongoing investigation by 
research groups in diverse disciplines ranging from Physicsfl, 2], to Geology [3, 4], and to 
Chemical Engineering [5, 6]. 

There are at least two, and possibly more, mechanisms thought to be at work: first, 
the swimmer can expand its hindquarters and thrust its wedge-shaped head into the bed 
[7], and second, it can submerge in an undulatory motion [8]. Neither of these behaviors 
is well understood, and both appear to represent active manipulation of the properties of 
the bed to make it solid enough to provide traction, but fluid enough to minimize friction - 
indeed, recent investigations of the motion of robots on sand beds hinge on understanding 
this balance between solid-like and fluid-like properties of granular beds [9] . 

To start with a simple model, we have performed simulations based on a recently proposed 
"Pushme-Pullyou" mechanism [10]. As shown in snapshots (Fig. 2, right), this mechanism 
consists of contracting and protruding the forequarters while enlarging the hindquarters 
(snapshots at t = 0, T), then enlarging the forequarters to anchor the new position and 
contracting the hindquarters (t = 2T, 3T), and finally returning to the initial state (t = 
4T, 0). We find that this mechanism does produce steady forward motion through a granular 




FIG. 1: A snapshot of sand swimmers {Plestiodon reynoldsi) emerging from a granular bed. Figure 
credit: Alessandro Catenazzi. 
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FIG. 2: Geometry of the swimmer (left) and its shapes at the end of each cycle of the 
swimming stroke (right). Movies of the swimming motion appear at http://bopper.t.u- 
t okyo .ac.jp/" shimada /sands wim /movies . 



bed. This mechanism differs from the well-known Purcell's swimmer which functions in 
viscous fluids [11-15], however it is amenable to direct simulation, provides insight into 
mechanisms that may be at work in granular systems [cf. SA Koehler et al., preprint 2008], 
and may enable intriguing future technologies such as reconnaissance or exploration for 
landmines beneath the desert floor. 

To simulate this system, we define a two dimensional swimmer consisting of anterior and 
posterior disks with diameters Da and Dp, connected by a linear spring with spring constant 
K and natural length L Q (t) (see Fig. 2, left). Both the swimmer disks and the surrounding 
granular particles have hard core potentials. The particle dynamics are modeled using an 
event-driven simulation method as described in [16-18]. To make a swimming stroke, we 
let the swimmer disks inflate and shrink, and correspondingly we define the natural length 
of the connection between neighboring disks to extend and retract. We consider reciprocal 
swimming strokes consisting of four cycles (A, B, C, and D, as illustrated in Fig. 2 right) 
with equal period (T) as: 
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where t' = (t mod AT). The initial shape of the swimmer is given as 
(L (0), D A (0), D P (0)) = (Lq, D-, D + ). AD = D + — D~ is the difference between the max- 
imum diameter D + and the minimum diameter D~ of the swimmer disks. AL = Lq — Lq 
is the difference between the maximum length Lq and the minimum length Lq. We set the 
bond to be very rigid (the spring constant K ~ 10 4 ) so that the actual length L(t) tightly 
follows Lo(t) in the sand media. In our work here, we fix the proportions of the swimmer to 
be Lq /Lq = D + / D~ = 2 and Lq = 2.5D~, hence the swimming motion is specified by the 
maximum diameter D + and the stroke time AT. Note that the scaled swimming frequency 
F = D + /AT gives the characteristic velocity of the swimming motion. We set the density of 
the swimmer at P sw j mmer = 0.68, so that the swimmer has approximately neutral density 
with respect to the surrounding sand with volume fraction 0.68 (loosely packed case such 
as rapid swimming under free-boundary) to 0.82 (system with compacting top wall) and its 
particle density p sanc i = 1.0. 

The sand particles are polydisperse, defined to have diameters uniformly distributed on 
1.0 ± 0.15 and we set the restitution coefficient to be e ganc j = 0.7. Collisions between sand 
particles and the swimmer are perfectly elastic ( e sw i mme r = 1-0)- To avoid inelastic collapse 
[19], we apply the TC model [20]: if a particle has a collision in a very small time window 
after its last collision, the collision is treated as elastic collision (e ganc j = 1). In this study 
we set t c = 10 -3 . We normalize the simulations by setting gravity, g, to unity, and we define 
the bottom of the system to be a fixed hard wall with specular reflection. The left and the 
right boundaries are periodic. For the top, we have studied two different conditions: a free 
boundary and a wall with constant pressure. In the latter case, we put a rigid and horizontal 
wall with specular reflection on the top with line density p wa ]\ = 100. We use this condition 
to mimic swimming in the "bulk", i.e. far beneath the surface. 

When we have a free-boundary on the top, dips are formed on the sand surface by the 
Pushme-Pullyou swimming strokes (Fig. 3). In this situation, the swimming direction is 
unstable because of the interaction between the swimmer and the free surface shape. The 
voids produced by the swimmer modify the free surface and change the stress field around 
the swimmer in a complex way. Therefore we focus on swimming in the more straightforward 
"bulk" condition. It is worth mentioning that we can modulate the swimming direction by 
adding quick deflation and moderate inflation motion. The surrounding sand cannot catch 
up with fast motion therefore the disk drops to the bottom of the void shaped by the disk 
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FIG. 3: A typical snapshot of the system with free boundary on top. Note that as the swimmer 
moves, it produces void spaces that propagate upward and modify the free surface (see text). 

itself before the deflation. In contrast, the disk almost keeps its position during the slower 
inflation cycle. Hence the disk get downward shift in total. However, a better control of 
direction is beyond the scope of this study. 

Providing a heavy top surface stabilizes the swimming direction. So we can evaluate 
swimming speeds and efficiencies as a function of the model parameters. In Fig. 4, we show 
the time evolution of velocity and power expended by the disks with moderate swimming 
frequency (F = 0.125), in which only small voids are observed around the swimmer. From 
the time evolutions shown, the swimming motion can be understood as follows. The small 
velocity of the larger disk in cycle A (posterior disk) and cycle C (anterior disk) implies 
that the larger disk is effectively anchored and moves little as the smaller disk extends and 
retracts. We have confirmed that the drag force is independent of velocity by simulations 
at multiple different frequencies. This is characteristic for slow motion of large intruders in 
granular media [21]. 

In contrast, during the inflating cycle (cycle B for the anterior disk and cycle D for the 
posterior disk) the inflating disk is held at a constant pressure p, i.e., the sand around the 
swimmer appears to become fluidized and consequently the power that the inflating disk 
exerts grows proportionally to time. This picture changes during the shrinking cycle at the 
same frequencies: the relaxation of the surrounding sand particles cannot catch up with the 
swimmer's motion, hence the swimmer cannot retrieve the energy it has expended. From 
this observation, one can roughly estimate the total displacement Ax and the total work W 
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FIG. 4: Time evolutions of the power and the velocity of the disks in the bulk condition. During 
cycle A, the smaller anterior disk is extended; during cycle B, the anterior disk expands and 
the posterior one contracts, during cycle C the posterior disk is retracted, and during cycle D 
the anterior disk contracts and the posterior one expands. The parameters for this example are 
D+ = 20 and T = 40 (F = 0.125). 

done by the swimmer during one cycle to obey: 

Ax ~ AL oc D + , (2) 
W ~ ^{( J D + ) 2 -( J D-) 2 }+2r /j D-AL oc^, (3) 

where p and 77 are phenomenological pressure and dynamic drag coefficient, and S = 
7r(.D + ) 2 /4 is the maximum area of the swimming disk. 

Let us now focus on two important quantities of the motion: the swimming velocity 
V = Ax/4T and the swimming efficiency E = SAx/W. The velocity and the efficiency as 
functions of swimmer size and frequency are shown in Fig. 5. In the slow-motion region, 
the swimming velocity increases linearly with the swimming frequency F since the total 
displacement obtained by each stroke is constant. As the swimming frequency is further 
increased, the resulting swimming velocity reaches a maximum and finally decreases. The 
decrease in swimming speed coincides with an enlargement in void space surrounding the 
swimmer: it appears that when the agitation is sufficiently rapid, a low density bubble forms 
around the swimmer, preventing it from gaining traction and moving forward. As a result, 
there is an optimum frequency for swimming velocity. 

What is more striking is that there is a second, different, swimming frequency at which 
the swimming efficiency is maximal (Fig. 6). The swimming efficiency takes its maximum 
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FIG. 5: Plot of the swimming velocities. The line "no slip" corresponds to the velocity expected 
for the swimming without slip: V = 1.25F. 
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FIG. 6: Plot of the swimming efficiencies. All swimmers have optimum swimming frequencies. 



in the slow-swimming regime, in which the swimmer does not slip against the bed. This 
implies that the efficiency maximum comes from the minimum in the total work W. In Fig. 
7, we confirm this: the total work shows a peak near the point at which the peak of the 
swimming velocity appears (scaled frequency F ~ 2), after which the stroke starts slipping. 
In the low frequency limit, we expect the total work expended over a swimming cycle to 
approach a finite value, unlike the case of swimming in fluid, where the work approaches 
zero. Our simulation results indeed suggest this non-zero quasi-static limit. In addition, it 
is found that the total work W takes a minimum at certain low, non-zero frequencies. How 
and why does the work exhibit a minimum before the quasi-static limit? Looking at the 
work done by one disk in each cycle separately, we find that the dragging cycle (cycle A for 
the anterior disk and cycle C for the posterior disk) and the inflating cycle (cycle B for the 
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FIG. 7: Total work done by the swimmer with diameter D+ = 20 scaled by S (bold line) and the 
corresponding scaled works done in each cycle by the posterior disk (thin lines). Note that the 
posterior disk anchors, deflates, is dragged, and inflates in cycle A, B, C, and D, respectively. 

anterior disk and cycle D for the posterior disk) contribute most to the minimum (Fig. 7). 
In contrast, the work done by the anchoring disk (cycle C and A for the anterior and the 
posterior disks, respectively) and the work retrieved to the deflating disk (cycle D and B for 
the anterior and the posterior disks, respectively) converge to a static limit around F — 0.1. 



This can be understood in more detail by scrutiny of the power expended by the swimmer. 
As displayed in Fig. 4 and argued by Eq. (3), we find the linearly increasing form in the 
inflating cycle under conditions that the disk feels constant pressure from the surrounding 
sands. As frequency decreases, the profile shifts to a rectangular shape which is similar to 
the dragging cycle. The profile in the dragging cycle keeps its rectangular shape for lower 
frequency, while increasing its height. Since the temporal evolution of the position of the wall 
on top is independent of the swimming frequency in the slowly swimming region, the global 
density of the sand is considered to be constant. Thus, the increase of the drag experienced 
by the disk as it deforms the surrounding sand comes from the re-organization of force chains 
within the granular bed. In other words, our simulations indicate that sand swimming relies 
on simultaneous solidification (near the anchored disk) and fluidization (near the moving 
disk) of the granular bed. This in turn relies on the swimmer's ability to exploit multiple 
timescales by moving one disk rapidly while the bed remains slightly fluidized, and holding 
the other disk static in its solidified environment. The swimmer at the efficiency optimum 
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takes advantage of this difference in timescales to minimize its energy expenditure. We 
speculate that more realistic models of sand swimming may depend on a similar balance 
between timescales in the granular bed. 

In conclusion, we have studied swimming in sand under Avron's Pushme-Pullyou mech- 
anism [10] using an event-driven simulation. We find that both swimming velocity and 
swimming efficiency have optima at different swimming frequencies. If the swimmer moves 
too fast, the voids created by the swimming stroke become larger and finally surround the 
whole swimmer. In this situation the swimmer slips and therefore the swimming velocity 
exhibits its maximum in fast motion region. If the swimmer moves too slowly, the bed 
re-solidifies before the swimmer can move forward. Optimal efficiency is achieved when the 
swimmer achieves a balance between these extremes: the inflating and moving disk is sur- 
rounded by weakly fluidized grains, while the anchoring disk is surrounded by solidified sand. 
These results suggest that the interplay between coexisting fluidized, slightly fluidized, and 
static granular states may be important for the understanding more complex situations in- 
cluding swimming near a free surface, the interaction between several swimmers [15, 22, 23], 
and more realistic shapes of swimmers. The frequency and size of the swimmer, i.e. its 
characteristic parameters, determine the time-scale over which the fluidization-solidification 
transition can be induced. In that sense our swimmer constitutes a simple device to locally 
probe this fundamental time-scale in a granular packing. 
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